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Abstract 

In this work we propose a generalization of Winfree's abstract Tile Assembly Model (aTAM) in which 
tile types are assigned rigid shapes, or geometries, along each tile face. We examine the number of distinct 
tile types needed to assemble shapes within this model, the temperature required for efficient assembly, 
and the problem of designing compact geometric faces to meet given compatibility specifications. Our 
results show a dramatic decrease in the number of tile types needed to assemble n x n squares to 0(yiog n) 
at temperature 1 for the most simple model which meets a lower bound from Kolmogorov complexity, 
and O (log log n) in a model in which tile aggregates must move together through obstacle free paths 
within the plane. This stands in contrast to the 0(log n/ loglog n) tile types at temperature 2 needed in 
the basic aTAM. We also provide a general method for simulating a large and computationally universal 
class of temperature 2 aTAM systems with geometric tiles at temperature 1. Finally, we consider the 
problem of computing a set of compact geometric faces for a tile system to implement a given set of 
compatibility specifications. We show a number of bounds on the complexity of geometry size needed 
for various classes of compatibility specifications, many of which we directly apply to our tile assembly 
results to achieve non-trivial reductions in geometry size. 
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1 Introduction 



The stunning diversity of biological tissues and structures found in nature, including examples such as 
signaling axons stretching from neurons, powerfully contracting muscle tissue, and specifically tailored coats 
protecting viral payloads, are composed of basic molecular building blocks called proteins. These proteins, 
in turn, are assembled from an amazingly small set of only around 20 amino acids. So how is it that so much 
structural and functional variety can be derived from so few unique components? The simplified answer 
is "geometry" . Essentially, a protein's function is determined by its 3-dimensional shape, or geometry. 
The exact sequence of amino acids which compose a protein (along with environmental influences such as 
temperature and pH levels) determine how that particular string of amino acids will fold into a protein's 
characteristic 3-dimensional structure. However, as simple as it may sound, the resulting geometries are 
often extremely complex, and predicting them has proven to be computationally intractable. It is from such 
geometrically intricate structure that nearly all of the complexity of life as we know it arises. 

Scientists and inventors have always recognized nature as providing invaluable examples and inspiration, 
and as for many other fields, this is also true for the study of artificial self-assembling systems. Self-assembling 
systems are systems in which sets of relatively simple components begin in disconnected and disorganized 
initial states, and then spontaneously and autonomously combine to form more complex structures. Self- 
assembling systems are pervasive in nature, and their power for creating intricate structures at even the 
nano-scale have inspired researchers to design artificial systems which self-assemble. One such productive 
line of research has followed from the introduction of the Tile Assembly Model (TAM) by Winfree in 27 . 
As a basic model, the TAM has proven powerful, providing a basis for laboratory implementations [5|[7|l6 
19)[20j[23)[24','29' as well as copious amounts of theoretical work [6l|9)[ll|[l2)[l8j[^|26j[28] . However, in this 
work, we have once again looked to the guidance provided by nature, this time in terms of the power and 
importance of the geometric complexity of the components of self-assembling systems, to extend the TAM 
in an attempt to harness that power. 

1.1 Overview 





Figure 1: The use of jigsaw faced macro tiles for self-assembly is emerging in both theoretical and experi- 
mental work. This figure contains three separate recent examples. The first figure depicts the theoretical 
technique of encoding binary strings within the geometry of tile growth into the third dimension, as seen by 
small blue tiles in this figure [9]. The second figure depicts a macro tile assembled using staged assembly 
from smaller tile types [10[ . Finally, the third figure depicts the experimental work of [14' in which a jigsaw 
geometry on the face of tiles is created with the DNA origami technique. 
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Table 1: Summary of our Results, tr denotes tlie number of distinct glues of a tile system to be simulated, with (7n and cr^ 
denoting only the number of north/south and west/east glue types respectively. L{M) is the size of the smallest geometry that 
can satisfy a binary n X m compatibility matrix M. 



We introduce a generalization of the abstract Tile Assembly Model (aTAM) in which tile types are 
assigned rigid shapes, or geometries, along each tile face. This model is motivated by the plausibility of 



implementing novel sophisticated nanoscale shapes with technology such as DNA origami 21 . We show 
that this model permits substantially greater efficiency in terms of tile type complexity when compared to 
assembling shapes in the basic temperature 2 aTAM. Furthermore, these efficiency improvements hold even 
at temperature 1. 

1.2 Results 

The abstract tile assembly model (aTAM) 27 , as well as many of the nanoscale self-assembly models 
spawned by it, feature single stranded DNA sequences as the primary mechanism for decision making. This 
commonality applies to weak systems such as deterministic temperature- 1 assembly, as well as stronger ones 
that rely on higher temperatures or stochastic methods. Since it is known that DNA strands are capable 
of hybridizing with sequences other than their exact Watson-Crick compliments, it is therefore reasonable 
to consider a tile assembly model in which one glue can potentially bond with an arbitrary subset of the 
other glues, with possibly differing strengths. Aggarwal et. al. have shown that such a non-diagonal glue 
function allows for significant efficiency gains in terms of the numbers of unique tiles used to assemble a target 
shape. Despite this potentially promising result, it is also true that designing non-specific hybridization pairs, 
while possible, is severely limited in a practical sense, and would likely introduce a potential for error in a 
much greater sense than is already present in laboratory experiments. 

If non-specific binding is impractical or impossible to implement, but powerful in theory, the question 
remains: are there any other mechanisms by which this power can be realized? One possible answer to 



this question is motivated by advances in DNA origami 14 21 in which DNA strands can be folded into 
blocks with semi-rigid jig-saw faces (see the rightmost image in Figure [T]). In this work we introduce a 
generalization of the aTAM in which tile faces are given some rigid shape (which we hereon refer to as 
geometry) . As suggested in Figure [2] the geometric hindrance which can be provided by this geometry 
is capable of simulating non-diagonal glue functions by creating a set of compatible and non-compatible 
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Figure 2: Examples of geometric tiles. Note that only the black portions on the corners are binding surfaces 
with glues, while the "teeth" in between provide potential geometric hindrance. Left: Compatible tiles. 
Right: Incompatible tiles (colliding teeth, which prevent the glue pads from coming together, are circled). 



faces. We show that this new model realizes much of the power of non-specific hybridization. Among our 
results, we show that n x n squares can be assembled in Q{y/]ogn) distinct tile types, which meets an 
information theoretic lower bound for the model and improves what is possible without geometric tiles from 
Q(\ogn/ log log n) (see [22]). In addition, this tile efficient construction requires only a temperature threshold 
of 1, thus showing this model can mimic both non-specific glue functions and temperature 2 self-assembly 
simultaneously. 

Next, we show that temperature-1 systems utilizing geometry can efhciently simulate a powerful class 
of temperature-2 aTAM systems. This class of systems, called zig-zag systems, is capable of simulating 
arbitrary Turing machines and therefore universal computation. Furthermore, the simulation performed 
using geometric tiles is efficient in that it requires no increase in tile complexity (i.e. the number of unique tile 
types required) or in the size of the assembly. This is especially notable due to the fact that it is conjectured 
(although currently unproven) that temperature-1 systems in the aTAM are not computationally universal 
(see [T3|[T5] for more discussion about temperature-1 assembly in the aTAM). 

While tile geometries provide a method for greatly reducing the tile complexity required to build squares 
in a seeded model similar to the aTAM (i.e. one in which tiles can only combine with a growing assembly 
one at a time), our next result holds for geometric tiles considered within the 2-handed assembly model 



(sometimes referred to by other names) [l|[2j|4j[8 [10 18[|28| . We show that, in this model, the tile complexity 
required to build a square is reduced to only O(loglogn) tile types. The construction presented utilizes 
the ability of 2-handed assembly to grow assemblies by the combination of sub-assemblies composed of 
groups of previously combined tiles, and, coupled with complex geometric patterns on the tile edges, forces 
assembling components to undergo intricate patterns of relative motion in order to combine with each other. 
The tile geometries required are, however, complex (O(lognloglogn)) and in a 2-dimensional model require 
disconnected components. We then show a simple extension to 3 dimensions which allows for connected 
components while retaining all other features. 

Finally, we conduct a detailed analysis of problems related to computing necessary patterns for tile 
geometries given specifications of the desired compatibility matrices (i.e. the listings of which tile sides 
should be compatible and incompatible with each other), with the goal being to minimize the size of the 
necessary geometries (as well as the running time of the computations). They deal with designing tile face 
geometries as a subset of Zi x Zi. Their solutions help show the feasibility and limitations of geometric 
tile face designs. We show a number of lower and upper bounds related to variants of the problem, some of 
which are incorporated into the previously mentioned constructions. 

1.3 Organization of this paper 

The remainder of this paper is organized as follows. In Section [2] we describe and define the new models 
introduced here. In Section [S] we present our constructions and proofs related to the self-assembly of n x n 
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squares using 0(-\/log n) tile types, as well as the simulation of zig-zag, temperature-2 aTAM systems by 
temperature-1 systems with geometric tiles. Section [4] describes our construction which utilizes geometric 
tiles as well as 2-handed assembly to self-assemble n x n squares using O (log log n) tile types. Additionally, 
there is a technical appendix which contains the majority of the proofs and construction details for the 
results presented in the previous sections, as well as the results related to computing compatibility matrices. 



2 Model 

In this section we define the basic geometric tile assembly model (GTAM) and the two-handed planar geo- 
metric tile assembly model (2GAM). We begin with an informal description of the aTAM. We then define 
the Geometric Tile Assembly Model (GTAM). The GTAM generalizes the aTAM ^ by adding a geometry 
to each tile face that may prevent two tiles from attaching. 



2.1 Basics 

A tile type is a unit square with four sides, each having a glue consisting of a label (a finite string) and 
strength (0, 1, or 2). We assume a finite set T of tile types, but an infinite number of copies of each tile type, 
each copy referred to as a tile. A supertile (a.k.a., assembly) is a positioning of tiles on the integer lattice 
Z^. Two adjacent tiles in a supertile interact if the glues on their abutting sides are equal. Each supertile 
induces a binding graph, a grid graph whose vertices are tiles, with an edge between two tiles if they interact. 
The supertile is r-stable if every cut of its binding graph has strength at least r, where the weight of an 
edge is the strength of the glue it represents. That is, the supertile is stable if at least energy r is required 
to separate the supertile into two parts. A seeded tile assembly system (TAS) is a triple T — (T, r, s), where 
T is a finite tile set, r is the temperature, usually 1 or 2, and s € T is a special tile type denoted as the seed. 
Given a TAS T — (T, r, s), a supertile is producible if either it is the seed tile, or it is the r-stable result of 
attaching a single tile r € T to a producible supertile. A supertile a is terminal if for every tile type r € T, 
r cannot be r-stably attached to a. A TAS is directed (a.k.a., deterministic or confluent) if it has only one 
terminal, producible supertile. Given a connected shape X C Z^, a TAS T produces X uniquely if every 
producible, terminal supertile places tiles only on positions in X (appropriately translated if necessary). 



2.2 Geometric Tiles and the Basic Geometric Tile Assembly Model (GTAM) 

In this paper we generalize the basic aTAM by assigning a geometric pattern to each side of a tile type along 
with its glue. For each tile set in the GTAM, fix two values w,£ £N. While at a high-level we still consider 
tiles as occupying unit squares within the plane, in order to determine whether or not adjacent tiles are 
geometrically compatible with each other, we define a tile body to be an x £ square (see Figure [s]), and we 
define a (tile face) geometry to be a subset of Z„ x Z^. A geometric tile type consists of a tile body which 
has both a glue and a geometry assigned to each side. For a tile type t, let northGeometry(t) denote the 
geometry assigned to the north side of i. Define eastGeometry{t) , southGeometry(t) , and westGeometry{t) 
analogously. Intuitively, the geometry of a tile type face represents the positions of inflexible bumps, or 
"filled-in" locations of the w x £ rectangle, that can prevent two tiles from lining up adjacently to one 
another so that the rectangles of their adjacent geometries completely overlap. Only if the w x £ geometries 
on adjacent sides of two combining tiles can completely overlap so that no location contains a filled-in portion 
of both, can any glues on those adjacent sides interact. Formally, we say a tile type t is east incompatible with 
tile type r if eastGeometry{t) P| westGeometry^r) ^ 0. We define north, south, and west incompatibility 
analogously. Seeded Geometric Tile Assembly takes place in the same manner as in the aTAM, with the 
added requirement that a tile type cannot be attached to a supertile at a position in which the tile type is 
either east, west, north, or south incompatible with another adjacent tile type in the supertile at a position 
west, east, south, or north, respectively, of the attachment position. As in the original aTAM, tiles are not 
allowed to rotate and must always maintain their pre-specified orientation, even while moving into position 
to attach to an assembly. 
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2.3 Two-Handed Geometric Tile Assembly Model 



The Two-Handed Geometric Tile Assembly Model (2GAM) extends 
the GTAM by allowing large assembled supertiles to attach to one 
another. We further restrict the model to planar assembly in which 
two supertiles may only attach if there exists a collision free path for 
the supertiles to traverse to reach their point of connection. In two di- 
mensional assembly this enforces that supertiles must be able to slide 
into position while staying in the 2D plane. With standard aTAM 
tiles, this requirement enforces that individual tiles of a supertile do 
not collide with individual tiles from another supertile while the su- 
pertiles shift into position. With geometric tiles, we must also enforce 
that the geometries of individual tiles do not overlap with other tile 
geometries. 



Informal Definition of the 2GAM As in the GTAM, tiles are 
composed of tile bodies and tile face geometries as shown in Figure [3j 
Within the 2GAM, two tiles may attach if 1) there exists a collision 
free path within the 2D plane to shift the tiles into an adjacent position in which the east (or south) geometry 
box of one tile exactly overlaps the west (or north) geometry box of the second tile, and 2) the east (north) 
and west (south) glues of each tile are equal and have strength at least r. More generally, preassembled 
multiple tile supertiles may come together if there is a collision free path in which the supertiles line up to 
create a r-stable assembly. The set of producible supertiles within the 2GAM is defined recursively: As a base 
case, all singleton supertiles consisting of a single tile are producible. Recursively, for any two producible 
supertiles a and /3 such that there exists a collision free path within the plane to shift a and /3 into a r-stable 
configuration 7, then the supertile 7 is also producible. The subset of producible assemblies of a 2GAM 
system to which no producible assembly can attach defines the terminally produced supertiles. Intuitively, 
this set represents the set of assemblies we expect to see from a system if it is given enough time to assemble, 
and we refer to this as the output of the system. A 2GAM is directed (e.g., deterministic, confluent) if 
it has only one terminal, producible supertile. Given a connected shape AT C Z^, a 2GAM F produces X 
uniquely if every producible, terminal supertile places tiles only on positions in X (appropriately translated 
if necessary). 

Please refer to Section|^for a more formal definition of the 2GAM model. Additionally, for a discussion of 
the different types of tile face geometries that are possible and the classes into which they can be categorized, 
please see Section [B] 
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Figure 3: Definition of a geometric tile. 



3 Complexities for the GTAM: Squares and r = 1 Assembly 

In this section we examine the power of the GTAM in the context of efficiently building squares and simulating 
temperature t — 2 ATAM systems at r = 1. We first show in Secton [O] that the tile complexity of n x n 
squares in the GTAM is Q{y'\og n) for almost all n by providing an order \/\ogn tile complexity upper bound 
construction for all n, and a matching information theoretic lower bound for almost all n. In addition, our 
upper bound construction utilizes only temperature r = 1. This stands in contrast to the temperature r = 2, 
0(logri/loglog?i) tile complexity result that can be achieved in the ATAM |3 . 

As the square construction shows, the GTAM seems to be powerful at t = 1. In Section 3.2 we consider 
the problem of simulating t = 2 ATAM systems within the GTAM, but at r = 1. We show that for a large 
class of temperature r = 2 ATAM systems called zig-zag systems, such a simulation is possible with no scale 
up in tile complexity or assembly size. Of particular note is the fact that zig-zag systems are capable of 
simulating universial Turing machines, something that is conjectured to not be possible in the ATAM at 
r = 1. 
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3.1 The Tile Complexity of GTAM squares: e{y/\ogn) 



In this section we analyze the size of the smallest tile type GTAM system that uniquely assembles an n x n 
square. Our first result is a construction that will assemble an n x n square using 0{\/\og n) tile types. We 
then show that this is tight for almost all n by applying an information theoretic argument to show that 
for almost all n, at least ^[s/Xogn) distinct tile types are required to uniquely assemble an n x n square. 
This result stands in contrast to the 8 (log n/ log log n) tile complexity for building squares in the standard 
ATAM model, showing that the GTAM is strictly more powerful than the ATAM. Further, our upper bound 
construction uses only temperature 1, while the ATAM construction requires temperature 2. 
In the remainder of this section we prove the following theorems: 

Theorem 3.1. The minimum tile complexity required to assemble an nxn square in the GTAM is 0{^\og n). 
Further, this complexity can be achieved by a temperature t — 1 system with 0{\/\ogn) size geometry. 

Theorem 3.2. For almost all integers n, the minimum tile complexity required to assemble an n x n square 
in the GTAM is n{^/\ogn). 

For the sake of brevity, we only give a high level overview of the upper bound construction and place the 
details in referenced appendix sections. 



3.1.1 0{\/\ogn) Construction Overview 

The tile system for the assembly of n x n squares in the GTAM at temperature r = 1 and tile complexity 
0{^\ogn) starts with the assembly of a roughly 2 x log n rectangle (2 x [log(^^)] + 2 to be precise) that is 
used as a base to encode a roughly logn digit binary number ( [log(^^iyi)] digits to be precise). The rectangle 
is efficiently built with 0{^J\ogn) tile types by using a tile set for simulating a 2-digit, base-ylogn counter. 
Such counters are known to exist in the ATAM at r = 2. To achieve r = 1 in the GTAM, we apply the 



transformation described in Theorem 3.6 to convert a zig-zag (see Definition 3.4 1 version of the r = 2 counter 
into a T = 1 GTAM system. The tileset for the basic r = 2 ATAM counter is given in Figure [TOj and the 
T — 1 GTAM version is given in Figure [TT] The details of this portion of the construction are described in 
Section [CMl 

The next step of the construction grows a third row of tiles on top of the surface of the roughly 2 x log n 
rectangle. Within this assembled row one tile type representing a binary or 1 bit is placed at each position, 
thus assembling a length roughly logn binary string upon completion. To generate O(logn) bits from only 
0{y/\ogn) distinct tile types is impossible in general within the ATAM. Within the GTAM, at this stage in 
the assembly we make use of the non-specific hindrance property of geometric tile faces to select the correct 
bit (and reject the wrong bit) at each of the logn bit positions. The key idea is that a collection of m 
geometric tile faces can be designed such that each face is compatible with a specified subset of the other 
faces, while incompatible with all others. Thus m tile faces can encode a compressed binary pieces of 
information (compatible or not compatible), thereby providing the possibility for a more tile type succinct 
assembly of an n x n square. A description of the decoder tile set is given in Figure [4] along with an example 
of how the assignment of geometry to tile faces permits the decoder tiles to efficiently select the correct bits. 
The details of this portion of the construction are described in Section C.0.3| 



Once the binary string is assembled on the surface of the 3 x log n rectangle, we utilize a well known 



constant sized set of tiles that implement a binary counter in the ATAM at t = 2 22 . This system reads 
a given surface of glues that denote an initial binary value for the counter, and then assembles upwards, 
incrementing a binary value encoded in tile types at every other row of the assembly. Once the counter is 
maxed out the construction stops, thus growing a rectangle of height roughly 2'°s" minus the initial value 



of the counter. This t — 2 ATAM counter construction is a zig-zag construction (see Definition 3.4). Thus, 
our construction applies Theorem |3.6|to convert to a r = 1 GTAM version. 



Finally, with the ability to generate large length 0(n) rectangles with 0{^/log n) tile types at t = 1, we 
combine 3 of these constructions to assemble the border of an n x n square using a factor of 3 times more 
tile types. A high level schematic of the approach is given in Figure |9] With the shell constructed, the 
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Figure 4: (a) The decoder set of tiles grows across the north face of the assembled counter from Figure |ll| to build a 
binary string of tiles by placing at each position a tile type representing either a 0-bit or a 1-bit. (b) For a given binary string 
B = 627 — 1 • • ■ b2bibo to be assembled, the geometries Gj, Gj, and Xi are assigned such that Gj and Xi are compatible if and 

only if b2r(7 i)—2ri+j — 0- i^) 3. given binary string B = b2r-l ■ • • b2bibo to be assembled, the listed compatibility matrix 

is obtained, along with the described assignment of geometry to each tile face. Such compatibility constraints can be achieved 
with geometry of length r for a length 2r^ — 2 string B. For more compact geometries in the case of less complex strings B, 
see Section [G| (d) An example assembly of decoder tiles. 



completion of the final rectangle can seed a growth of filler tiles to fill in the body of the square, finishing 
the construction. The final details of the construction are described in Section fC. 0.41 



3.1.2 Tight Kolmogorov Lower Bound for GTAM Squares: fl{y/\ogn) 

Theorem 3.3. For almost all integers n, the minimum tile complexity required to assemble annx n square 
in the GTAM is n{y/]ogn). 

Proof. The Kolmogorov complexity of an integer n with respect to a universal Turing machine U is Ku{n) = 
min \p\ s.t. U{p) = bn where 6„ is the binary representation of n. It is known that Kjj[n) > [logn] — A for 



at least 1 — {-^)^ of all n (see 17 for results on Kolmogorov complexity). Thus, for any e > 0, Kij{n) > 
(1 — e) logn = fi(logn) for almost all n. 

Consider a tile simulator program (of constant size in bits) that reads as input a GTAM tile system 
(encoded as a bit string). Suppose the simulator is modified so that it outputs the maximum extent (i.e. 
width or length) of a shape that is terminally produced by the input system. When such a simulator is 
paired with a GTAM system that uniquely assembles an n x n square, the combined program constitutes a 
program that outputs the integer n, implying that the total number of bits of the simulator (constant) plus 
the encoding of the tile set must be at least Kij{ri). As the simulator has a constant size, this implies that 
the number of bits to to encode the GTAM system must be at least Ku{n), which is r2(log n) for almost all n. 
To achieve our bound we now show that any GTAM system can be encoded using 0(|Tp) bits (independent 
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of the size/area of the tile face geometries) assuming a constant bounded temperature. To achieve this, we 
do not cxphcitly encode the geometry for each tile face, but instead utilize a compatibility matrix. 

Encoding a GTAM system. For a GTAM system F = (T, r, s), arbitrarily index each distinct face of 
each distinct tile in T from 1 to 4|T|. Define the compatibility matrix M for F to be the 4|T| x 4\T\ matrix 
such that Mij = 1 i is the index of an east (or north respectively) edge and j is the index of a west 

(south respectively) edge and i and j have incompatible edge geometries. M can be encoded using 0(T^) 
bits, and the remaining portions of F can easily be encoded in asymptotically fewer bits, yielding an O(T^) 
bit encoding for any GTAM system. Note that even without the explicit representation of the GTAM's 
geometries, a simulator can derive what the system will build from the compatibility matrix M. 

Now consider the smallest tile type GTAM system F = (T, r, s) that uniquely assembles an n x n square. 
As F can be encoded in 0(|rp) bits, we know that for almost all n, ci|Tp > C2logn for constants ci,C2- 
Therefore, |T| = il{\/\og n) for almost all n. □ 

3.2 Simulating Temperature r = 2 ATAM Systems with r = 1 GTAM Systems 
Definition 3.4. Zig-Zag System. A tile system F ~ {T, t, s) is called a zig-zag system if: 

1. The location and type of the i*^ tile to attach is the same for all assembly sequences. 

2. The i*'' tile attachment occurs to the north, west, or east (not south) of the previously placed tile 
attachment in all assembly sequences. 

For the proofs of the following two theorems and technical details about the notion of tile system "sim- 
ulation" , please see Section |D] 

Theorem 3.5. Any temperature r = 2 zig-zag ATAM tile system F = (T, 2, s) can be simulated by a t = 1 
GTAM tile system T = {R,l,q) with tile type scale |i?|/|T| = 0(1). The simulation utilizes geometry size 
at most loga„ -f- log log cr„ + 0(1) where (t„ is the number of distinct north/south glue types represented in 
T. 

Theorem 3.6. Any temperature r = 2 zig-zag ATAM tile system F = (T, 2, s) can be simulated by a r = 1 
GTAM tile system T = {R,l,q) using only 1 non-null glue type and tile type scale |_R|/|T| = 0(1). The 
geometry size of the simulation system is at most logu -I- log log cr + 0(1) where a is the number of distinct 
glue types represented in T. 



4 2 GAM Results 

In this section, we explore the theoretical limits achievable when utilizing geometric tiles by designing tiles 
whose edges contain highly complex geometries. Furthermore, we move to the 2-handed variant of the 
GTAM, the 2GAM, to allow for the geometric hindrances experienced by individual tiles to be grouped and 
combined to provide more complex interactions between larger supertiles. The goal, rather than providing 
a realistic and potentially experimentally realizable set of constructions, is to gain further understanding 
into the interplay between geometry and the types of computations which can be carried out via algorithmic 
self-assembly. 

We now present the details of our construction, which reduces the tile complexity required to self-assemble 
an n X n square to a mere O(loglogn) tile types, while requiring a geometry size of O(lognloglogn). Our 
construction requires the constraint of planarity, in which components are not allowed to float into position 
from above or below the assembly, but must always be able to slide into position with a series of translations 
along only the x and y axes. However, the intricate geometric designs and complex series of movements 



require that individual tile geometries are composed of disconnected components. (Note that in Section E.5 
we show how to extend the tiles into the third dimension, utilizing a total of 4 planes, in a manner which 
results in connected tiles and also implicitly enforces the restriction that only tile translations along the x 
and y axes must be sufficient to allow for tile attachments.) 



8 



4.1 Self-assembly of an n x n square with O(loglogn) tile types 

Theorem 4.1. For every n £ N, there exists a 2GAM tile systera F = {T,2) which uniquely produces an 
n X n square, where |T| — O(loglogn), and with O (log n log log n) size geometry. 

To prove Theorem |4.1[ we present the following construction. 



4.1.1 High-level sketch of the construction 



T 

n"+4 




Following is a list of values based on the partic- 
ular dimensions of the square to be formed and 
which are used throughout the following discus- 
sion: 

• n: dimensions of the square to self- 
assemble 

• n': [logn] 

• n": [logn'] 

• s: 2"' + 2"" + 2n" + 8-n 

• h: 2""-i 



1 



Figure 5: A high-level sketch of the construction for building a 
square in the 2GAM. 



• Ci : 2-handed counter which counts from 
through 2" — 1 for a total of 2" columns 

• C2: standard counter which counts from 
s through 2"' - 1 for a total of 2"' - s 
columns 

• C3: 2-handed counter with "buffer" 
columns which counts from through 
2"" - 1 for a total of 2""+i - 1 columns 



Figure [S] shows a high level view of the main 
components of this construction. Without loss 
of generality, we can consider the construction to be composed of a series of sub-assemblies, or modules, 
which assemble in sequence, with each module completely assembling before the next begins. The careful 
design of all modules ensures that none can grow so that they occupy space required by another, and that 
each will be able to terminally grow to precisely defined dimensions that result in the final combination 
forming exactly an n x 71 square. For the rest of this discussion, we will describe the formation of the 



modules in such a sequence. (See Figure 12 for a series of high-level images which exemplify the ordering of 
the formation of the square from these modules.) 

this construction makes use of one counter, 



Similar to the construction in Section 3.1 



Ci, to assemble 

an encoding of a number which in turn seeds another counter, C2. Ci assembles in a 2-handed manner, 
meaning that each number which is counted is represented by exactly one one-tile-wide column of tiles, 
and individual columns form separately and then combine to form the full counter of length 2" (similar in 
design to counters found in |12 ). Each column of the counter, besides representing a counter value, is used 
to represent (on the north face of the northernmost tile) one bit of the seed value s for C2 . Each column 
can form in one of two versions: one that represents a 0, and one that represents a 1, for the corresponding 
bit of s. The east and west sides of the tiles forming the columns of this counter contain geometries which 
force the columns, in order to combine, to "wiggle" up and down in patterns based on the counter values 
of those columns. See Figure [6] for an example pair of compatible columns. The columns also contain tiles 
with geometries which "read" those patterns of wiggling and allow columns to combine with each other if 
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and only if they are the correct versions of the counter columns, namely those with the seed bit values which 
correctly correspond to their location in the counter. It is the tiles of this component as well as those of 
the counter C3 to which the intricate geometries are applied, and thus they receive a much more detailed 
explanation in Section [E.l[ 

C2 is a standard binary counter (i.e. one that would also 
assemble correctly in the aTAM) which utilizes 16 tile types (see 



Figure 15 ) and grows to complete the majority of the western 
side of the square. Next, a small set of 7 "filler" tile types (see 
Figure 14b I fill in the majority of the square, and once they 



have filled in a sufficient portion of the northern portion they 
provide a platform to which C3 can attach (as long as C3 is fully 
formed). In order to provide a directed system with only one 
terminal assembly, the "incorrect" columns (those which couldn't 
become part of Ci) are able to combine into the 2-handed counter 
structure C3 via some extra buffer columns (see Figures 14 and 
19 1. Finally, the filler tiles are able to complete the formation of 
the square. Note that the tile types which make up C2 and the 
filler tiles require no geometries but only standard glues. 

By utilizing the assembly of supertiles (i.e. sub-assemblies of 
grouped tiles) and carefully designing geometries which force the 
supertiles forming Ci to move in well-defined patterns as they 
attach, we are able to essentially "transmit" information about 
tiles in one location of a supertile to the interfaces where poten- 
tial binding is occurring with other tiles in the same supertile. 
By concatenating this information from such a group of distant 
tiles, the binding "decision" can be made based on an arbitrarily 
large amount of information (as long as the geometry sizes scale 
appropriately). This results in a dramatic lowering of the tile 
complexity required to assemble an n x n square, with the trade- 
off being an increase in the complexity of the tiles themselves. 

Please see Section |E] for much more detail and several sup- 
plementary images describing this construction. Additionally, 
a comprehensive example has been provided in Section |E.4| to 
which the reader can refer for additional clarity. 

4.2 Analysis of tile complexity and geometry 
size 



l-Cap 
Tile 3 



1-Cap 
Tile 2 



1-Cap 
Tile 1 



1-Cap 
TileO 



0-Cap 
Tile 3 



0-Cap 
Tile 2 



0-Cap 
Tile 1 



0-Cap 
TileO 



First, we analyze the tile complexity of this construction, module 
by module, in order to determine the overall complexity. 
The tile complexity of each component is as follows: 

Counter tiles: There are n" bit positions which each 
require a constant number of tile types (as can be seen 
from the depiction in Figure 14a), plus the requirement for 



Fi gure 6: Example columns for the counter Ci. 
Note that all colored areas are fiUed-in, and areas 
colored white are empty, although they may be 
outlined for reference. 



a hard-coded column on each of the west and east sides, for a total of O(loglogn) tile types. 

Cap tiles: There are 4 cap tile positions which each 
need to be able to represent a or a 1 cap, for a total of 8 
tile types. 



C2: 16 tile types. 
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Counter tiles: There are n" bit positions in the buffer 
columns which each require a constant number of tile types, 
for a total of O (log log n) tile types. 

Buffer cap tiles: There are 4 cap tile positions which 
each require a single tile type, for a total of 8 tile types. 

• Filler tiles: 7 tile types. 

Thus, the total tile type complexity is 0(loglogn)+0(l)+0(l)+0(loglogn)+0(l)+0(l) = O(loglogn). 

Next, we simply note that the geometries defined for all tiles in this construction consist of rectangles of 
dimensions (2"" +h + 4)x (n" + 2) = (2"" + [2""/2] + 4) x (n" + 2) = 0(logn x log log n), and therefore 
the geometry size is O(lognloglogn). 
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A Formal 2GAM Definition 



Hindrance Map To permit a supertile to interweave itself into an attachable position requires modeling 
translations of a supertile at the resolution of the size of individual units of the tile face geometries. To this 
end we define a Hindrance Map for a supertile a that represents the set of positions that a supertile takes 
up, including the bodies of each tile in the supertile, along with the positions blocked by each geometry for 
each tile face. 

Formally, consider a supertile a. The Hindrance Map Ha is the following set of positions: For each tiled 
position {x^y) in supertile a, the following points are defined to be in Ha- (Please refer to Figure [t] for 
depictions of tile components in terms of w and ^, and note that the entire area occupied by a 

BODY: {{i,j)\x ■ [w + t) + w < i < x ■ {w + t) + w + I, 
y ■ (w + e) + w < j < y ■ {w + £) + w + 1} 

WEST GEOMETRY: {w + £) < i < x ■ {w + £) + w, 

y ■ {w + £) + w < j < y ■ (w + £) + w + I, 

{i — (x ■ (w + i)),j — {y ■ {w + £) + w)) G W estGeometry{a{x , y))} 

SOUTH GEOMETRY: {{i,j)\x ■ [w + £) + w < i < x ■ {w + + w + I, 
y ■ (w + £) < j < y ■ (w + £) + w, 

(i — {x • {w + £) + w), j — {y ■ {w + £))) e SouthGeometry(a{x, y))} 

NORTH GEOMETRY: {{i,j)\x ■ {w + £) + w < i < x ■ [w + £) + w + I, 
+ 1) • (w + ^) < j < (?/ + 1) • (w + £) + w, 
{i — {x ■ {w + £) + w),j — ((y + 1) • {w + £))) G SouthGeometry{a{x, y))} 

EAST GEOMETRY: {{i,j)\{x + l) ■ (w + £) < i < (x + 1) ■ {w + £) + w, 
y ■ {w + £) + w < j < y ■ {w + £) + w + I, 

(i — {{x + 1) ■ {w + £)),j — {y ■ {w + + w)) e WestGeometry{a{x, y))} 

The final hindrance map Ha is the union of the sets BODY, SOUTH, NORTH, EAST, and WEST for 
each tiled position of a. 

Planar Translation of Supertiles Given two supertiles a and (5, a collision free translation of /3 with re- 
spect to a is any translation of Ha that can be obtained by a sequence of unit translations {t^o, Wi, U2, W3, . . . Uj,} 
where vq is an initial translation that shifts Ha such that all positions of Ha are northwest of all positions 
of Hj3, and each m is one of the translations {u„ = (0, l),Ue — (—1,0), Ms — (0, — l),Utu = (—1,0)}. Further, 
after each translation Ui, it must be the case that Ha does not overlap Hp. A grid locked collision free 
translation is a collision free translation in which Ha has been shifted by multiples of w + £ in both the x 
and y direction. We are interested in grid locked translations as they correspond to direct translations of 
a at the resolution of tiles, rather than the higher resolution translations of Ha- We require grid locked 
translations for a supertile to attach to another supertile as such a translation is needed for the tiles to line 
up. 

2GAM Model A Two-Handed Planar Geometric Tile Assembly System consists of a duple (T, r) where T 
is a set of geometric tile types and t is the positive integer temperature of the system. Given a 2GAM system 
r — {T,t), a supertile is producible if either it is a single tile from T, or it is the r-stable result of a grid 
locked collision free translation of two producible assemblies. A supertile a is terminal if for every producible 
supertile f3, a and /3 cannot be r-stably attached. A 2GAM is directed (e.g., deterministic, confiuent) if 
it has only one terminal, producible supertile. Given a connected shape AT C Z^, a 2GAM P produces X 
uniquely if every producible, terminal supertile places tiles only on positions in X (appropriately translated 
if necessary). 
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Figure 7: Mapping of tile coordinates (red) to dimensions and units of tile bodies and geometries. 



B Classes of Tile Face Geometries 

Here we discuss different classes into which tile face geometries can be classified. Note that all GTAM 
results presented in this paper have "bump" geometries, while the 2GAM result (in its 2-dimensional form) 
has "unrestricted" geometries. See Figure [8] for an example of each class. 

Unrestricted This is the most general class of geometries and places no restrictions on 
the portions of a tile face geometry region (i.e. the w x £ rectangle) which are filled-in 
and which are empty. Such geometries may be infeasibl to implement as the pieces of a 
geometric face may not be connected to the tile body. However, given a third dimension 
it is plausible that such a scheme might be implemented by attaching particles to the face 
of a substrate which also attaches to the tile body. 



Connected A slightly more restricted class, the connected class allows arbitrary pat- 
terns to be filled-in in within the tile face geometries as long as all such portions retain a 
connected path to the tile body. 



Bar A bar geometry is restricted to lines of filled-in points which are connected to the 
tile body and extend directly away from it. Each bar can be of length x where < x < w. 



Bump Bump geometries are the simplest possible types of geometry and consist of a 
set of points which are directly connected to the tile body. This class can be thought of 
as simplified bar geometry with w — 1. 



Figure 8: Classes 
of tile geometries. 
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C Additional Details for the O^Iogn GTAM Square 
Construction 



C.O.l Construction Notation 

Consider a positive integer n (the widtli of tlie square we wish to assemble) . For the sake 
of clarity, assume n is even. 

• Let n' = [log^]. 

• Let r = r^(l/2)n' + l/2] . 

• Let B ^ 2"' - n/2 - 1. 



CO. 2 Base 0(-^log n) Counter at Temperature t = 1 

The first step of the construction is the assembly of a 2 x log n rectangle that will serve as a bed for a third 
layer that will place a row of tiles that represent a logn bit binary number. The r = 1 GTAM system used 
for this portion of the assembly is obtained by applying the transformation from Theorem |3.6| to convert an 
efficient r = 2 ATAM system into a equivalent r = 1 GTAM system. 

The temperature r = 2 ATAM system that will be converted is described in Figure [TO] and constitutes a 
2-digit, base r counter. The counter works within the ATAM at temperature r = 2 and is a generalization 



of the base-2 version first described in 22 . By specifying the east glues of tiles Si and S2, the counter can 
be seeded to any specified starting value. In the tile set given in the figure, the counter is initialized to value 
0. The value of the counter is incremented at every other column as the assembly grows from west to east, 
finally halting when the counter rolls over to 0. Thus, for a given choice of r and an initial seed value b, the 
final assembly will be a 2 x 2r^ — 2b rectangle. For our construction, we utilize r = f-^/ (l/2)ri' + 1/2] , which 
guarantees enough room to place a length n' binary string in the next step of the construction. As we can 
initialize the counter to be shorter if needed, we assume the counter has been initialized to grow to length 
exactly n' + 2 (an extra 2 positions are not used to encode bits in our constructions, thus the extra 2 length). 
To modify the ATAM system of Figure [TO] to a r = 1 system, we observe that it is a zig-zag system 



according to Definition 3.4 (rotated 90 degrees). Therefore, we can apply Theorem 3.6 to obtain an equivalent 
T ~ 1 GTAM system shown in FigurejlT] The general case details of the conversion are detailed in Section 3.2 
but the basic idea of the transformation is to replace east strength 1 glues with the null glue type, and assign a 
unique geometry to each edge for each glue type. In particular, for an east/west glue type x, a corresponding 
pair of geometries are computed, E(x) and W{x), such that E{x) is incompatible with all other geometries 
within the system with the exception of W{x)^ and vice versa. Each occurrence of the glue type x on the 
east face of a tile type is replaced by a GTAM tile type with geometry E{x) for the east face geometry. The 
same replacement by W{x) is done for west occurrences of x. This geometry assignment is also applied to 
all north/south glue types as well. The result is a system that assembles in the same fashion as the original 
temperature t — 2 system and with the same tile complexity, but does so at temperature t — 1. 



CO. 3 Bit Decoder Tiles 

The next step in the construction consists of a collection of decoder tiles which grow across the surface of 
the 2 X 2r^ rectangle assembled from the previous step. The general tile set for these decoder tiles is given 
in Figure |4] The growth of these tiles is initialized by the blue glue displayed by the final tile placement 
of the base r counter from the previous section. The decoder tiles consist of a repeating chain of 2r tiles 
with labels to 2r — 1. The west geometry of tiles with label i are compatible with the east geometries of 
tiles with label {i + 1) mod 2r, and incompatible with all other geometries. Thus, the chain of tiles must 
assemble in the proper order. Further, there exactly 2 tile types with each integer label from to 2r — 1, a 
green type and an orange type. Our goal is to assemble a supertile that encodes a given target binary string 
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Phase 3: Binary Counters to Build Shell of Square 
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Figure 9: This fi gure shows a high level overview of the different steps in the assembly of n X n squares in the GTAM with 
0{\/log n) tile types at temperature r = 1. In phase 1, a tile set that implements a 2-digit, base 0(\/log n) counter is used 
to form a length O(logn) bed upon which a binary number will be assembled. Phase 2 places green and orange decoder tiles 
which denote either a or 1 bit at each position of the third row of the assembly. In phase 3, the assembled binary string is 
utilized as the seed for a binary counter set of tile types which grow a length n rectangle. Phases 1-3 are repeated two more 
times to create the outer shell of an n X n square. Finally, a collection of filler tiles are seeded to fill in the body of the square. 
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Base r2-digit Counter at Temperature 2 
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Figure 10: This figure contains a tile system that assembles a 2 X rectangle for a given integer r. The construction is an 
implementation of a 2-digit, base-r counter with tile complexity 5r + 2 that assembles at temperature t = 2 in the standard 
ATAM. 
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Figure 11: This tile system is the result of applying the transformation from Theorem 3.6 to the zig-zag counter described 
in Figure [To| The systems utilizes a single (red) glue at temperature r = 1. For each east / west (or north/south) glue type x 
in the initial ATAM system, there exist two corresponding geometries W{x) and E(x) (or N{x) and S{x)) in the new GTAM 
system such that W{x) is compatible with E{x), but incompatible with all other geometries in the system. 
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along its surface. The encoding is the pattern of green and orange tiles with orange tiles representing binary 
bits, and green representing binary 1 bits. The key to get the goal binary string assembled is to enforce 
that at each bit position the correct bit is chosen. We do this by appropriately assigning geometry to the 
south face of the decoder tiles and the north face of the base r counter tiles. 

In more detail, suppose we are given a target binary string B = 62r-i ■ • ■ &2&i^o to be assembled, the 
geometries Gj, Gj, and Xi are assigned such that Gj and Xi are compatible if and only if b2r(j—i)-2ri+j = 0- 
A detailed example of such a compatibility matrix is given in Figure |4| along with a sample set of geometry 
assignments to tile faces that satisfies such constraints. More generally, such compatibility constraints can 



be achieved with geometry of length r for a length 2r — 2 string B by way of Theorem G.5 which is 



asymptotically the best achievable in most cases. In the case of less complex binary strings, more compact 
geometries can be obtained as discussed in Section [G] From these compatibility constraints, the desired 
target binary string is guaranteed to assemble. Further, the final placed tile can be specified as a special 
type with a north purple glue which seeds the next portion of the construction. 



CO. 4 Binary Counter Tiles 

The next portion of the construction, seeded by the final tile placed during the decoder tile portion, consists of 
a binary counter tile system which grows north, start from the binary string decoded in the previous section, 
up until the counter rolls over. The construction is a t = 1 GTAM version of a well known temperature t = 2 



ATAM construction of 0(1) tile types 22 . The counter increments every other row, and thus will grow to 
a height of exactly 2" +^ — 2B — 2, where B is the initial value of the counter which is the binary string 
encoded in the previous phase of the construction. Our goal is for the counter to build to a height equal to 
n, the dimension of the goal nxn square. We thus choose B in the previous section to be i? = 2" — n/2 — 1. 

Finally, the construction is finished by observing that the binary counter construction can be implemented 
such that the final tile placed at the northeast corner of the assembly exposes a glue type which seeds the 
assembly of a rectangle that grows east for exactly n — n' — 2 units. This can be accomplished in the exact 
same way we achieved a north growing rectangle of length n, but with an alternate choice of initial binary 
string assignment. This construction can in turn seed a south growing rectanle of the same length. This 
final rectangle can seed the growth of a final 0(1) size collection of tile types which fills in the casing of the 
hollow nxn square, yielding the final full nxn square. A high level figure depicting the construction is 
given in Figure |9] 



D Details of Zig-Zag Simulation 

Definition D.l. Nice Zig-Zag System. A zig-zag system F = (T, r, s) is called a nice zig-zag system if: 

1. The terminal assembly of F contains no exposed east-west non- glue types. 

2. All producible assemblies of F contain no mismatched glues. 

Definition D.2. Tile System Simulation. An ATAM or GTAM system F2 = (T2,r2,S2) is said to 
simulate a second ATAM or GTAM system Fi = {Ti,ti,si) if: 

1. There exists a function F : Ti — > P{T2) such that an assembly sequence ((ti, xi, yi), . . . (ti,Xi,yi)) is 
valid for system Fi if and only if there exists tile types g F{ti) such that the assembly sequence 
((ri G F{ti),xi,yi), . . . {n e F{U), Xi.yi)) is valid for F2. 

The tile complexity scale factor of the simulation is defined to be |r2|/|ri|. This definition only considers 
the case of simulating a system without scaling the size of the assembly (scale factor 1). See [9] for a more 
general definition of simulation that permits scaled assembly size factors. 

Observation D.3. The "simulate" relation between tile systems is transitive. Further, if system A simulates 
B with tile type scale factor x, and B simulates C with tile type scale factor y, then A simulates C with tile 
type scale factor xy. 
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Lemma D.4. Any zig-zag system Fi = (Ti, ti, si) can be simulated by a nice zig-zag system r2 = (T2, T2, S2) 
with tile type scale factor |T2|/|ri| = 0(1). 



Proof, of Theorem 3.5 Consider a zig-zag system T = (T, 2,s). By Lemma D.4 there exists a zig-zag 



system F' = (T', 2, s') that simulates F with 0(1) tile type scale such that the assembly of F' has the "nice" 
properties described in Definition |D.1| 

We now define a system T = {R, l,q) that simulates F' with 0(1) tile type scale, and thus simulates F 
with 0(1) tile type scale by the transitivity of simulation. 

Let a'^ denote the set of all west-east glues that are represented in the tile set T' . Let ctJj denote the set 
of all north-south glues that are represented in the tile set T' . Let H denote an injective mapping from the 
glue types represented in T' to some new set of strength- 1 glues p. 

For each t £ T' , define the the geometric tile type rt as follows. Denote the north, south, east, and west 
glues of t as north(t), south{t), east{t), and west{t) respectively. Let the west glue of rt be west(rt) = 
H{west{t)) and east{rt) = H{east{t)). If north{t) is a strength-2 glue, then north{rt) = H{north{t)), 
otherwise north(rt) = null. If south{t) is a strength-2 glue, then south{rt) — H{south(t)), otherwise 
south{rt) = null. 

We assign the empty set geometry to all east/west edges of tile types in R. We assign non-empty 
geometries to north/south edges of each tile type rt such that the south edge of a given rt G R is compatible 



with a north edge of a tile type r^ € R if and only if south(t) — north{v). By Theorem G.9 we can achieve 
such compatibility requirements with a geometry of size at most log cr^ -I- log log . 

We now show that the system T = {R, 1, rg') simulates F' = (T', 2, s') (with tile type scale factor 1). Con- 
sider the function F{t) = {rt} to map tile types from T' to R. Suppose T correctly simulates F' up to the first 
i—1 steps. That is, for the first i — 1 assembly sequence steps ((ti, Xi, yi), . . . (ii_i, x^-i, j/,;_i)) of the unique 
assembly sequence of F', the first i—1 steps of the assembly sequence of T are ((rt^ , xi, yi), . . . {rt-_-^, Xi^i,yi^i)) . 
Consider the z*'' step in the assembly sequence of F', {ti, Xi,yi). We know that the placement position {xi, yi) 
must occur north, west, or east of the previously placed tile position (x^-i, yi_i). If the placement occurs 
to the north, then (r^., Xi, yi) is a valid step in the assembly sequence for T. This is because the tile trans- 
formation explicitly assigns tile type a strength- 1 north glue that matches the south glue of r^. , as ti_i 
must have a north strength-2 glue that matches the south glue of ti. Further, by the the fact that F' is 
nice, there are no exposed east/west glue faces of the F' assembly after i — 1 steps, and thus the T assembly 
will have no other exposed glues (of strength-1) beyond the single north glue of rt. ^. This implies that the 
placement of rt- at position {xi,yi) is the only possible next tile placed for system T. 

Now suppose the i*'* tile attachment occurs to the east of the previously placed tile in the F' assembly. 
The required strength-2 attachment threshold for the i*'' tile can be achieved by two strength-1 glues, or by 
a single strength-2 glue. In the cooperative strength-1 glue case, we know that rt- is the only tile in R that 
both matches the east glue of rt-^-^ and has a compatible geometry with the north face of the tile type at 
position {xi, yi — 1). Thus, (rt^ , Xi, yi) is a valid next element of T's assembly sequence, and is the only valid 
next element by the nice properties of F'. For the case of a strength-2 east attachment, we know that the 
{rt-,Xi,yi) attachment is vahd for T because of the the matching strength-1 east glue of rt._-^ and west glue 
of rj. j , and because the tile type south of (xi, yi), if there is one, is guaranteed to have a compatible north 
geometry with rt- by the "no mismatched glue" property of nice zig-zag systems. The attachment is also 
unique because of the "no exposed glues" property of nice zig-zag systems. The remaining west attachment 
case is analogous to the east attachment case. 

Therefore, T = {R,l,rs') simulates F' = (T',2,s') with tile type scale |i?|/|T'| = 1, and therefore also 
simulates F = (T, 2,s) with tile type scale 0(1). The size of the geometry of the simulation is at most 



log cr^ -I- log log (T^ by Theorem G.9, and is thus log2cr„ -I- loglog2a-„ — logcr„ -I- loglogcr„ -I- 0(1) where cr,, 



is the number of north/south glue types in T. □ 
Proof, of Theorem |3.6| We use the same approach for simulating nice zig-zag systems as given in the proof 



for Theorem 3.5 with the modification that all non-null glue types from system T' are mapped to a single 
strength-1 glue type x. Further, rather than empty set geometry assigned to all east west glues in the 
simulation set, we assign geometries that satisfy a compatibility matrix which assigns a to entries which 
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correspond to identical glues, and 1 to entries corresponding to non-identical glues. Thus, while there is a 
single glue, only the appropriate tile with geometry representing the appropriate east /west glue will attach. 
By the same analysis given for Theorem |3.6[ we achieve a simulation set for any zig-zag system with 0(1) 
tile type scale and logo- -I- log log cr -I- 0(1) geometry size. □ 



E Additional details for 2GAM Results 

This section includes additional details and images describing the construction in Section |4] 




Figure 12: A sketch of the assembly sequence of the square and its components: 1. first the counter module Ci assembles; 
2. the initial value s is seeded by Ci, which allows C2 to assemble; 3. filler tiles complete a sufficient portion of the square; 4. 
A fully formed version of C3 attaches; and 5. filler tiles complete the square. 



E.l The 2-handed counters Ci and C3 

For ea ch til e in this construction, the bodies are squares of (2" +/i + 4)x(2" + 4) geometric units. (See 
Figure 13a for an example.) The north and south geometries consist of (2" + /i + 4) x {n" + 2) rectangles of 



geometric units (i.e. rectangles of length £ — 2^ + h + 4 and width w = n" + 2) . All north geometries are 
completely empty, while all south geometries are completely filled in. However, the east and west geometries 
are composed of {n" + 2) x (2" +h + 4:) rectangles of geometric units (i.e. rectangles of length £ = 2" +h + A 
and width w = n" + 2) which contain intricate collections of gaps, which we call sockets, and projections, 
which we call prongs. 

Note that while the techniques utilized in this construction can be generalized to form counters of arbitrary 
bases, the tile complexities are asymptotically identical regardless of the base, so for simplicity of explanation 
we utilize only base 2 counters. 

The purpose of the counter module Ci is to count from through 2" — 1, for a total of 2" values, each 
represented by exactly one column. The northern glue on the northernmost tile of each of column is used to 
represent one of the bits of the value s (from left to right, the most significant bit to the least significant). 
Each column is composed of n" + 4 tiles arranged vertically. The tiles in each column perform two tasks: 1. 
represent an n"-bit binary number m for < m < 2"' — 1, and 2. present a northern glue which represents 
the correct value (0 or 1) for the 2" — (rn + l)th bit of s. The southern n" tiles, which we call the counter 
tiles, each represent one bit of m, from top to bottom the most significant bit to the least significant. The 
northernmost 4 tiles are gadgets called caps which serve to represent the bit values on the north of the 
columns, while ensuring their correct positioning relative to the bits of s. There are two possible caps which 
can form, a 0-cap and a 1-cap. Either cap can nondeterministically attach to any column of counter tiles to 
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Figure 13: Details of the geometries of the counter tiles. 



provide the 4 northernmost tiles of any column (except for the leftmost and rightmost columns, which are 
special cases discussed later). This allows for the formation of two versions of the columns which represent 
each counter value m: one that represents a on the north and one that represents a 1. It is the purpose 
of the cap to ensure that only the version of the column with the correct bit of s (i.e. for the column 
representing the value m and thus the 2" — (m -I- l)th bit of s, the one representing the bit value matching 
the 2" — (m + l)th bit of s) can attach to neighboring columns. 

Each column can completely form independently of every other column, and in fact must be fully formed 
before it can combine with any other column. For each bit position of the counter, there are two tile types, 
one representing a and one a 1. Only those for the least significant bit position have strength- 1 glues on their 
east and west, while those for every other bit position have null glues on their east and west sides. However, 
all counter tiles have geometries on their east and west sides which will be explained shortly. The north and 
south sides of these tiles have flat geometries and strength-2 glues which allow them to nondeterministically 
combine with the bit values above and below (but only in the correct order of significance). Thus, the 
counter portions of each column can form by nondeterministic combinations of exactly one tile type for each 
bit position, and therefore columns can form which represent each of the values through 2" — 1. For each 
counter value, two columns can form: one with a 0-cap and one with a 1-cap. The east and west sides of the 
cap tiles have geometries but null glues, except for the northernmost which has strength-1 glues on the east 
and west. The north and south sides of each cap tile are similar to those of the counter tiles in that they 
have only flat geometries and strength-2 glues, except for the northernmost cap tile which has a strength-1 
glue that represents the bit of s provided by that column. (The rightmost, and therefore least significant, 
bit of s will be represented by a strength-2 glue, and both the most significant and least significant bits will 
be represented by special glues which convey both the bit values and the most and least significance of those 
positions.) The positioning of the two strength-1 east and west glues in each column, on only the top and 
bottom tiles, ensures that the only way that two columns can possibly combine with each other (keeping in 
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(a) Template for the tile set which allows the formation of columns representing binary 
numbers in the range 1 through 2^ — 2 {preventing the formation of numbers and 
2" — 1), and which increment the binary numbers represented by each column from 
west to east. The bit significance positions for each group of tiles are shown. Note 
that a unique group of tile types is created for every bit position, although those for 
positions 2 through n" — 2 arc shown in a single group. Each tile has strength-2 glues 
on its north and south edges (except for the north edge of the northernmost group 
and the south edge of the southernmost group). The east and west bit values are 
represented in the actual tile set by the east and west geometries. 



(b) Filler tiles which fill in the bulk 
of the square. 



mind that this construction operates at temperature 2) is if they are fully formed (i.e. they contain all 4 cap 
tiles and exactly one counter tile for each bit position). 

E.1.1 Counter tiles 

Figure [T4a| shows a tile set which performs portions of the functionality of the counter tile types. This tile set 
is designed to form 1 x n" tile columns representing rt"-bit binary numbers, with specific groups of tile types 
designed for various groups of bit positions. All tiles have strcngth-2 glues on their north and south sides 
(except for the north glue of those in the most significant bit position and the south glue of those in the least 
significant bit position) which match only those of the adjacent edges of tiles for adjacent bit positions. The 
columns formed represent n"-bit binary numbers on their west and east sides such that for each number x 
represented on the west, the value x + 1 is represented on the east. Additionally, every possible n"-bit binary 
value can form except for the values and 2" — 1 (i.e. all O's or all I's). Those values are represented by 
the leftmost and rightmost columns, respectively, of the counter Ci and are formed by a separate set of hard 
coded tile types since the leftmost column must expose a special glue on the north which marks that bit of 
s as the most significant bit, and the rightmost column must expose a special strength-2 glue on the north 
which marks that bit of s as the least significant bit, and it must also have flat geometries and strength- 1 
glues along its eastern side which can attach to the filler tiles. 

For each bit position p (where < p < rt"), in order to represent bit values bi on the west and 62 on 
the east, the geometries are designed as shown in Figure [T3a[ The 1x2 rectangles labeled bi and 62 each 
have exactly one position filled in and one left empty, depending on the bit values being represented and 
following the patterns shown in Figure |13b[ It is these regions which enforce the restriction that only if two 
columns represent the same bits in all positions can they line up and come together completely to combine. 
The remaining area of the west geometry is completely empty except for one unit located h + 3 units from 
the southern edge and p + 2 units from the east. This single point is referred to as the prong, which will 
slide into a complementary socket in the east geometry of a tile in the same bit position of a column to the 
left. 
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The portion of the east geometry to the left of region 62 is a 2 x rt" + 1 rectangle which is completely filled 
in. The geometry immediately to the north of that creates the socket, which has one empty column p + 1 
units from the east and one additional empty spot h units up that column and one unit to the west. The 
area labeled B2 consists of p + 1 columns which are divided horizontally into a number of blocks dependent 
upon p, the bit position being represented. It is divided into 2" ~p blocks of width p + 1 and height 2^. 
Those blocks can be thought of as being numbered from the southernmost as 0, through the northernmost 
as 2" ~P. If 62 ~ 0, then each even numbered block is empty and each odd numbered block is filled in. If 
62 = 1, then the even numbered blocks are empty and the odd numbered blocks are filled in. Examples can 
be seen in Figure [17] 

Intuitively, the eastern socket is split into two groups of regions - those representing a and those 
representing a 1. If the socket is part of a 0-tile (i.e. a tile which represents a for the particular bit of the 
counter value), the regions representing are left open and those representing 1 are filled in. For a 1-tile, 
the opposite is true. The number of regions is determined by the bit position. For the most significant bit 
position, there are only two regions, with the bottom representing and the top 1. The tile for the next bit 
down has four regions, which are essentially each of the regions of the previous tile split into 2 sections, with 
the bottom of each new pair representing and the top 1. This splitting and doubling continues downward 
until the tile for the least significant bit which has 2" regions consisting of individual squares. In a 0-tile, 
the bottom of each pair starting from the bottom is left open and the top of each is filled in, and vice versa 
for a 1-tile. 

The portion of the east geometry above the socket, which is an empty n" + 2x h rectangle, is referred to 
as the padding region. A padding region exists on the north of each west and east geometry, and its purpose 
is to guarantee that the geometries of tiles in one bit position of a column will not collide with the geometries 
of those in adjacent bit positions of a column correctly combining with it. This padding region is sufficient 
because the maximum relative translation that two columns which can combine with each other will ever 
need to experience is a vertical offset of 2""/2, and the padding region plus the separation provided by the 
north-south geometries of combined tiles is at least that large. 

E.1.2 Cap tiles 

The remainder of each column of Ci consists of four tiles called the cap which combine to form a vertical 
column one tile wide and 4 tiles high. There are two possible types of cap, a 0-cap and a 1-cap, which can 
nondeterministically attach to any partial column representing a counter value to form the top four tiles of 
a complete column. It is the job of a 6-cap (for 6 = or 1) to ensure that only if it has attached to a counter 
value y (which will be positioned at the location of the (2" — {y + l))th bit of s) and b matches that bit 
of s, only then can that column attach to other columns within Ci (specifically, those columns representing 
y — I and y -I- 1) and thus become a portion of that counter. Otherwise, an unmatching column (i.e. one 
to which the "wrong" cap has nondeterministically attached) is relegated to use within C3, which will be 
described later. 

Please refer to Figure [14] for a graphical depiction of the geometries of the cap tiles. The north and south 
geometries of all cap tiles are the same as those of the counter tiles: a completely filled in (2" -I- ft, -1-4) x {n"+2) 
south geometry and a completely empty (2" + h + A) x (n" -I- 2) north geometry. The prong in the west 
geometry of cap tile is located h + 5 units from the souther and the prong in east geometry of cap tile 1 is 
located 2h + 3 units from the north. Each consists of a single filled-in location 3 units away from the bodies. 
The remainder of those geometries are empty. Additionally, the west geometry of cap tile 2 and the east 
geometry of cap tile 3 are completely empty. Thus, it is the 4 sockets, one on each cap tile, which provide 
the bulk of the functionality of the caps. 

The east geometries of cap tiles and 2 each have two empty rows on the south and a buffer region on 
the north. In between are sockets which each have one completely filled-in row above and below them. The 
west geometries of cap tiles 1 and 3 each have one filled-in row on the south, a buffer region on the north, 
two empty rows south of the buffer, one filled-in row south of that, and sockets in between the two filled-in 
rows. The socket patterns of the cap tiles occupy a single column and are referred to as Si and S2, plus 
their complements Si and 6*2. (Note that the sockets Si and S2 are not "utilized" during the combination 
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of columns which form into Ci, but only for those which become part of C3.) These columns are 2" units 
tall, and by associating each of these units with one of the 2" possible values of a counter column to which 
a cap can attach, it is possible to utilize the pattern of gaps and blocks in a socket to determine which type 
of cap (0 or 1) can be attached to a particular counter value while allowing that column to become part of 
Ci. By virtue of the fact that the prongs are relatively short compared with those of the counter tiles, and 
due to the patterns of sockets and prongs on the counter tiles, compatible columns (which can only combine 
if they are completely formed, so we only consider completely formed columns in this discussion) are forced 
to align in such a way that, during their combination, the prong of a cap which is incident upon the socket 
of a neighboring cap will be at a position which corresponds to the value represented by the counter value. 
In this way, we can allow the prong to slide into an "accepting" gap if and only if the cap type is correct for 
that value, namely if that counter value should be associated with a or 1 cap and thus a or 1 bit of the 
number s. 

To provide this functionality, the socket patterns for a cap of type b are the following. 5*1 treats the 2" 
units as representing each value /e for < fc < 2" — 1 counting downward from the top to bottom. For the 
fcth unit, if k coincides with a numbered bit position of s which has the bit value b, there is a gap, otherwise 
it is filled in. Si is the exact complement (which is equivalent to designing Si for bit value 1 — 6). 6*2 treats 
the 2" units as representing one "invalid" value plus each value A; for < /c < 2" — 2, counting upward 
from the bottom to top with the first position being the invalid, and thus always filled-in, position. Again, 
for the fcth unit, if k coincides with a numbered bit position of s which has the bit value b, there is a gap, 
otherwise it is filled in, and again, S2 is the exact complement (with the exception that for both the invalid 
location is filled in). 

Intuitively, the reason for the north-south reversal and offset of 1 for socket values between the pairs of 
sockets 5*1 and S2 is that the prongs that will be validating the combination of counter-to-cap values will 
be positioned (during the combination of a pair of columns) relative to a single counter value, and thus the 
sockets need to use that one value to validate the two columns which represent (in the correct case) two 
consecutive values. Additionally, the motion of the prongs relative to the respective sockets will be vertically 
reversed due to the fact that they are on opposite sides of the combining columns. (Also, recall that whenever 
two columns attempt to combine which do not encode consecutive values, the patterns encoded in regions 
foi and 62 will prevent their combination.) Essentially, S2 and S2 are designed so that if they are included 
in a column which encodes the counter value y, the positioning of the prong from a complementary column 
to the right, which would represent the counter value y + 1, would align it with holes if and only if the bit 
value of the {{y + 1) — l)th, or yth, bit of s equals b. 

E.1.3 Buffer columns and C3 

As previously mentioned, every counter column can form in two versions, one with a 0-cap and one with 
a 1-cap. The design of the glues and geometries is sufficient to guarantee that a fully formed Ci can and 
will include exactly one version of each column, namely that which will present the correct value for the 
corresponding bit of s. In order to ensure that this construction is directed, and thus has exactly one terminal 
assembly, the versions of columns which are not included in the formation of Ci must be included somewhere 
within the eventual nx n square. That is the purpose of C3 . 

Intuitively, C3 is a counter of approximately double the length of Ci which incorporates all of the "bad" 
columns (meaning those with the wrong caps to be included in Ci), but since they cannot directly combine 
with each other, each pair is separated by a buffer column that not only allows them to combine but ensures 
that they are in fact bad columns (and hence the approximate doubling of the length as compared to Ci). 

Each buffer column is formed of two portions, similar to the other columns. However, the bottom portion 
of buffer columns simply represent a single binary value from the range 1 through 2" — 1, inclusive. This 
value is represented on both the west and east sides rather than an incremented value appearing on the ciast 
side as in the counter columns, and the geometries which represent each bit are the same as those for the 
counter tiles. The logic provided by the north-south glues must simply ensure that columns representing 
every binary number from 1 through 2" — 1 can form, thus simply excluding a column of all O's. (This is 
because the counter column representing encodes the number 1 on its east, and thus the buffer column 
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need never combine with a column representing the number 0.) Thus, there will be a buffer column which 
can be sandwiched by every consecutive pair of bad columns. 

There is only one version of each of tlie buffer cap tiles, which can be seen in Figure [141 Similar to the 
tiles of the counter columns, the only non glues on the west and east edges are strength-1 glues on the east 
and west of the buffer counter tiles in the least significant bit position (the southernmost in each column) 
and the buffer cap tile 3 (the northernmost in each column). As shown, the only "interesting" geometry is 
that of the prongs on buffer cap tile 2 and 3. These prongs are positioned so that they must fit into the 
sockets of complementary counter columns which allow combinations of counter values with cap values that 
are the opposite of those required for inclusion in Ci. In this way, full buffer columns can combine on both 
sides with exactly those full counter columns which received the "wrong" cap and in the correct order to 
count from through 2" — 1 and therefore construct C3. 

E.2 Details of the counter C2 

Once Ci has completely formed, the binary number s is encoded on its northern face as a series of strength-1 
glues, along with a strength-2 glue on the north of the east-most tile, which is the least significant bit of s. 



To these glues, tiles of the types shown in Figure 15 attach to perform the counting of s through 2" — 1 



This is a basic counter which can assemble in both the standard aTAM as well as the 2HAM, resulting in 
the same produced assembly since it is polyomino-safe (see [28]). It is notable that as each row assembles, 
it "checks" whether or not it is the maximal value of the counter and, if so, causes a tile to be placed in the 
west-most position which has a special north-facing glue that will eventually allow C3 to attach. 

E.3 Completion of the square 

Figure [12] shows the sequence in which the full assembly can progress. Once Ci has completely formed and 
C2 has completely grown northward from it, filler tiles are able to complete a sufhcient portion of the rest 
of the square that a fully formed (and only a fully formed) C3 can attach. The filler tiles are shown in 
Figure [Mb] Once C3 has attached, the remaining filler tiles are able to complete the formation of the full 
n X n square. 

E.4 An example: constructing a 250 x 250 square 

Here we provide an example of portions of this construction as applied to a 250 x 250 square. 

• n: 250 

• n': [logn] = [log 250] = 8 

• n": [logn'l = [log 8] = 3 

• s: 2"' -I- 2"" + 2n" + 8-n^2^ + 2^ + 2(3) + 8 - 250 = 28 = OOOIIIOO2 

• h: 2""-i - 1 = 22 - 1 = 3 

• Ci : 2-handed counter which counts from through 2" — 1 = 2^^ — 1 = 7 for a total of 2^^ = 8 columns 

• C2: standard counter which counts from s = 28 through 2" — 1 = 2^ — 1 = 255 for a total of 
2"' - s = 256 - 28 = 228 columns 

• C3: 2-handed counter with "buffer" columns which counts from through 2^ — 1 = 7 for a total of 
2""+i - 1 = 2"^ - 1 = 15 columns 

Figure [16] shows the dimensions for the components of the square for this example, as well as the dimen- 
sions of the tile bodies and geometries of counter tiles. 
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Fi gur6 14: Dimensions of the features in the geometries of the counter column cap tiles as well as the buffer column cap 
tiles. White portions indicate areas which are not filled in, while grey areas are filled in (although socket areas Si, S2, Si, and 
S2 contain a mixture 26 
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Figure 15: The tile types which form the counter C'2, grouped by significance of bit position (i.e. tile types for the most 
significant bit position are in the leftmost group, those for the "interior" bit positions are in the middle group, and those for 
the least significant bit position are in the rightmost group. 
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(a) Dimensions of the modules of the square. (b) Dimensions of the geometries of the counter tiles. 

Fi gure I61 Dimensions of the components of a 250 x 250 square in the 2GAM construction. 
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Table 2: Patterns for each combination of socket and cap type. Each "X" represents a location which is 
filled-in. Note that the numbering for the Si sockets starts at the north and increases southward, while it is 
the opposite for the S2 sockets. 



The geometries of the counter tiles are a straightforward application of the generic definitions for the 
counter tile dimensions, which are based on the value of n" and h, to form versions which represent and 1 
for each bit position. However, the geometries of the cap tiles are also dependent upon the binary value of s 
(padded to 8 bits), in this case "00011100", since they perform a mapping of counter values to bit values of 
s. In this case, the geometries of socket 5*1 in the 0-cap, which has positions for all 8 possible counter values, 
must have openings in exactly the positions 0, 1,5,6, and 7 (and hence filled-in locations at positions 2,3, 
and 4). Recall that the positions are numbered from north to south for the Si sockets. The 1-cap should 
have the exact opposite pattern of openings and filled-in locations for its socket , as should the socket Si 
of the 0-cap. Clearly, the socket Si of the 1-cap should match the 0-cap's socket. 

The geometries of the S2 sockets perform a mapping of counter values to bit values of s as well, but they 
perform a "one-off" mapping by matching bit values of s with the counter values which are 1 greater than 
the bit positions. Additionally, the direction in which the positions of the socket are labeled is reversed, 
meaning they are counted from the south to north. Therefore, the sockets 6*2 for the 0-cap and 5*2 for the 
1-cap should have openings at exactly positions 1, 2, 3, and 7 counting from top to bottom, while the sockets 
S2 for the 0-cap and 5*2 for the 1-cap should have openings at exactly positions 4, 5, and 6. Note that the 0th 
position of ^2 sockets, due to the one-off mapping, is an invalid position which we always fill in (intuitively, 
because the counter value used to align with an 5*1 socket is always 1 greater than the number represented 
by the column containing the socket). Please see Table [2] for a full listing of the socket patterns and Figures 
[iTjjTSj and [T9| to view examples of counter and buffer columns. 



E.5 Equivalence of disconnected 2-dimensional tiles and 3-dimensional tiles 
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Figure 20: Side view of equivalent but con- 
nected 3-dimensional geometric tiles. 



The construction in this section has been presented as a 2- 
dimensional construction meaning that the tiles are 2-dimensional 
objects, and due to the constraint of planarity they are never allowed 
to travel out of the plane during assembly. However, the geometries 
defined for the tiles often contain disconnected portions, namely 
filled-in units which are not connected via an unbroken path of addi- 
tional filled-in units to their tile bodies. It is notable that by merely 
relaxing the restriction of forcing tiles to remain 2-dimensional, and 
instead allowing them to extend into three additional planes, this 
construction can achieve connectivity while removing the explicit 
restriction that translations of tiles during assembly must only oc- 



cur within the x — y plane since that restriction instead becomes implicit due to the design. See Figure 20 
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Figure 17: Two counter columns combining during the self-assembly of a 250 X 250 square (part 1). a) Two fully formed and 
correct counter columns representing the values 5 and 6, b) Translation of the right column to allow prong-to-socket alignment 
for the most significant bit (highlighted), and c) Translation for alignment of the second bit (highlighted). 
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Fi gure 18; Two counter columns combining during the self-assembly of a 250 X 250 square (part 2). a) After alignment of the 
final bit, all 4 prongs (highlighted) are able to slide into the back column of the sockets in the left column, b) Final north-south 
translation aligns each prong with the west-most gaps in each socket (highlighted) and also aligns the bit patterns of regions 
bi and 62 of each pair of counter tiles (also highlighted), c) After the final east-west translation, the counter columns are fully 
combined, allowing the glues to interact and stably bind them. 
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Fi gure 19l Columns formed during the self-assembly of a 250 X 250 square, a) A counter column with the wrong cap, which 
prevents its combination with the next consecutive counter column, b. The incorrect counter column can combine with a buffer 
column encoding the correct counter value. 
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for a depiction of how the tiles are extended along the z-axis. For each of the east and west geometries a 
solid filled-in square of units is added in the plane above and below, respectively. These attach to every 
unit of the geometries as well as to additional filled-in squares above and below the tile body to ensure 
full connectivity. Finally, an additional layer is added which covers the tile body and lies above the west 
geometry, with exactly one open plane between it and the west geometry. In this way, it is ensured that all 
portions of the tiles and their geometries are connected, and exactly the tiles that could connect with each 
other in the 2-dimensional planar version can now connect with each other. 



F Two Function Problem 

In this section, we study the problem of bar-to-bump reduction for the design of the tile face geometries, 
namely converting bar geometries to bump geometries (see Section [b]). The goal is to attempt to simplify 
the geometry patterns as much as possible. Called the two function problem^ it takes an input integer 
and requires the design of two functions / : {0, 1, . . . , ^} {0, 1}" and g : {0, 1, . . . ,1} — {0, 1}" to satisfy 
f{x) AND g{y) — 0" if and only if a; -I- y < /, where AND takes as arguments two binary strings, x and y, 
of matching length and returns the string of binary values representing the logical AND operation between 
each corresponding pair of bits in x and y. The target is to make n as small as possible. It is trivial to 
design the two functions with n = I. We show a lower bound that exactly matches the upper bound. 

Definition F.l. Let l,ne N. Define / : {0, 1, . . . , -> {0, 1}" and 5 : {0, 1, ...,;} ^ {0, 1}" where / and g 
take a number between and inclusive, and return a binary string of length n such that f{x) AND g(y) = 
0" if and only ii x + y < I. (Note that we define AND as the function which takes as arguments two binary 
strings, x and y, of matching length and returns the string of binary values representing the logical AND 
operation between each corresponding pair of bits in x and y. In this case, there must not be a 1 in the same 
position of both bit strings.) Goal: For a given I, find / and g such that n is minimal. 



We have Theorem F.2 for a lower bound for the two function problem. This lower bound exactly matches 



the upper bound shown in Theorem F.3 



Theorem F.2. For each I, every solution for the two function problem needs n > I. 

Proof. Assume that we have a solution for the two function problem. Let / : {0, 1, . . . ,1} — {0, 1}" and 
g : {0,l,...,l} ~> {0, 1}" such that f{x) AND g{y) = 0" if and only iix + y<l. 

For a string s in {0, 1}", define C(s) to be the number of Is in s. For example, C(OlOOll) = 3. Define 
D{f ) = Yl\=o C'(/(i)). For each string s — aia2 . . . a„ in {0, 1}", define G(s) = {i : Ui ~ 1 and 1 < i < n}. 
For each string s — 0102 . . . a„ in {0, 1}", define s[i] = a^. 

Let the solution / and g satisfy that D{f) is the largest for the least n. We claim that G(/(0)) C 
G(/(l))cG(/(2))c...cG(/(0). 

Assume that there is a m € {0,1, ... ,1} with G{f{m)) % G{f{m + 1)). Let i be the index such that 
f{m)[i] = 1 and f{m + = 0. Let f(m) = aia2 . . . a„ and f{m -I- 1) = &162 ■ ■ • b„. Define function /' as 
follows: f'{k) = f{k) for each fc 7^ m+ 1, and f'{m+ 1) = hi . . . bi^iOibi^i . . . 6„ = 61 . . . 6i_il6i+i . . . 6„. We 
claim that /' and g form a new solution for the two function problem. 

For each y with {m + 1) + y < I, we have m + y < I. Therefore, /(m) AND g{y) — 0". This implies 
o-i ■ 5(2/) W = 0- Furthermore, we also have /(m -I- 1) AND g{y) — 0". Therefore, 

C{f'[m + 1) AND g{y)) < Cif{m + l)ANDgiy)) + a,-g{y)[i] 

< + 
= 0. 

Thus, G(/'(m + 1) AND g{y)) = 0. Therefore, 

/'(m + 1) AND ff(y)=0". 



32 



For each y with (m + 1) + y > ^, we have /(m + 1) AND g{y) ^ 0". Since — 1, we have 

Cifim + 1) AND giy)) > Cif{m + l)ANDgiy)) 

> 1. 

Therefore, 

/'(to + 1) ANDg{y)^0". 

Combining the last two cases, we have that /' and g form a new solution for the two function problem. We 
have D{f) < D{f'). This contradicts the fact that D{f) is the largest for the same n. 

Therefore, the solution / and g has that G(/(0)) C G(/(l)) C G(/(2)) C . . . C G{f{l)). It is easy to see 
f{i) =^ f{j) ioTO<i < j <l because f{i) AND g{l - z) = 0" and /(j) AND g{l - i) ^ 0". Thus, we have 
that G{f{i)) is a proper subset of G{f{i + 1)). Therefore, n> I. □ 

Theorem F.3. The two function problem always has a solution with n = I. 

Proof. Let f{i) = I'O'"' and g{l - i) = On'"*. It is easy to verify that f{x) AND g{y) = O' if and only if 
x + y<l. □ 



G Matrix Problem 

In this section, we study the complexity of a matrix problem that is related to tile geometry design for 
compatibility specifications. An efficient solution for this problem can achieve reductions in geometry sizes. 

Called the matrix problem, the goal is generating a binary matrix via vector inner products. Given a 
0—1 matrix M of size m x n, find a list of vectors Ei, . . . , Em and a list vectors Wi , ■ ■ ■ , Wn of length 
I for each vector such that M{i,j) = if and only if Ei ■ Wj — 0, where • is the inner product of two 
vectors. The target is to minimize the length I, which is denoted by L{M) ~ I for the least I. We show 
that for almost all n x n binary matrices M, L{M) = <d{n). A submatrix Mi of M is characterized by 
{Ri,Ci) where Ri is a subset of row indices of M and Ci is a subset of column indices of M. In the 
case that all 1 entries are in submatrices Mi — Ci), • • • , Mf = {Rt, Ct) that satisfy Ri \^ Rj = % and 
Ci O Cj = ^ for i 7^ j, we show the equality L(M) — L{Mi) + • • • + L{Mt). This relationship gives a tool 
for solving a class of concrete matrix problems. For another case that all entries of M are in submatrices 
Ml = d), • • • ,Mt = {Rt,Ct) that satisfy Ri n i?^ = and d n Cj = for i ^ j,we derive the lower 
bound and upper bound max(L(Mi), • • • , L{Mt)) < L{M) < (1 + e) max(L(Mi), • • • , L{Mt)) + O(^) for 
an arbitrary constant e > using a randomized algorithm. This randomized algorithm can be used to find 
a matrix design with a small length The dimension of each submatrix Mi is vrii x in table [T] 

For two lists . . . , Em and Wi, . . . , Wn such that each Ei or Wi are of length I, they are called a l-list 
pair. 

Matrix Problem: Given a m x n — 1 matrix M, find a list vectors Ei, . . . , Em and a list vectors 
Wi, . . . , Wn of length I for each vector such that M{i,j) = if and only if Ei ■ Wj — 0, where • is the inner 

product of two vectors such that Ei ■ Wj = X]i=i ^kWk for Ei — 6162 • • • e; and Wj ~ W1W2 ■ ■ ■ wi. Minimize 
the length /. 

For an integer matrix Almxn, its rank is the number of linear independent rows. A matrix Mnxn is 
diagonal one matrix if all the diagonal elements are one and all other elements are zero. A matrix Mnxn is 
diagonal zero matrix is all diagonal elements are zero and all other elements are one. 

G.l Some Lower Bounds 

In this section, we show some results for the lower bound of the matrix problem. Most of the lower bounds 
results are for concrete matrix problems. 
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Definition G.l. For a binary matrix M, define L{M) to be tlie least lengtli I for a solution of the matrix 

problem M. 

Definition G.2. Assume that M = {aij)mxn is a binary matrix. 

• Define S{M) = min(m, n). 

• A submatrix M' — (R, C) of M is characterized by a set R of rows and a set C of columns in M . The 
elements of M' are all elements aij with i G R and j € C. We also define S(M') = min(|_R|, |C|). 

• Two submatrices Mi = (i?i,Ci) and M2 = {R2,C2) of M are independent if -Ri fl i?2 = and 

Ci n C2 = 0. 

• A set H of submatrices of M is independent if every two of them are independent. 

• Let H = {Ml, . . . , Mt} be a set of independent submatrices, define T{H) = max(S'(Mi), . . . , S{Mt)). 

• For a sequence s = aia2 ■ ■ ■ a„, and a subset P = {ii, . . . , i^} with ii < Z2 < . . . < of integers of 

{1,2, . . . ,n}, define s[P] = Ui-^Ui^ . . -Cii^. For an integer i < n, define s[i] = a^. 

Theorem G.3. Let Mi, . . . , Mf be independent submatrices of Mmxn and all Is of M are in Mi, . . . , Mj. 
Then L{M) = L{Mi) + ... + L{Mt). 

Proof. Assume that Ei, . . . and Wi,...,W„ form a solution of least length I for the matrix M = 
{ai j)mxn- Let Mi = {Ri,Ci) for i = ,t. Define Pj = {j : for some {u,v) G Ri x Ci with Eu[j] = 

WS] = !}• 

Claim 1. Eri [Pi], . ■ ■ , E^^ [Pj] and [Pj], . . . , Wc^ [Pj] form a solution for Mi, where Ri = {ri, . . . , r^} 
and Ci = {ci,...,c^}. 

Proof. By the definition of Pi, for each {u,v) E Ri x Ci, we have Eu[x] = or = for any v ^ Pi 

(otherwise, x is in P,). Thus, E^ AND Wy = O' if and only if E^[Pi\ AND Wy[Pi\ = O'-, where k = \Pi\. □ 

Claim 2. Pj n P,- = for i 7^ j. 

Proof. This can be proved by contradiction. Assume PiC\Pj 7^ 0. Let x G PiHPj. There is a fi) G P.; x 
such that Eui[x] = Wy^[x] = 1. There is a {u2,V2) G Rj x Cj such that £'M2[a;] = W^2[a;] = 1. Therefore, 
{Eu^ AND Wy^) contains a bit Eu-^ [x] ■ W^^ [x] = 1. Since Mi,. . . ,Mt are independent, ui G Ri and V2 G Cj, 
we have Oiui,v2 

= 0. This contradicts that {Eu, AND W^J = if and only if au,,v2 =0. □ 

Claim 3. L{M) > L{Mi) + ... + L{Mt). 

Proof. By Claim 1, since i?i[Pi], . . . , Em[Pi] and Wi[Pi], . . . , W„[Pi] form a solution for the matrix problem 
M, for i = 1, • • • ,i, we have L{Mi) < \P,\. By Claim 2, we have |Pi| + . . . + \Pt\ < I = L{M). Thus, 
L(M) = \Pi\ + ... + \Pt\> L{Mi) + ...+ L{Mt). □ 

Assume that Ei^n, ■ ■ ■ , Ei^ry and Wi^a, • • • , W^i.c^ form a solution for Mi = {Ri, Ci) with Ri = {n, . . . ,ry} 
and Cj = {ci, . . . , c^} for i = 1, . . . , The length of each vector for Mi is li = L{Mi). 

For J = 1, . . . , m, define = ^1 . . . At such that Aj = O'' if z ^ PiU. . . Pj, and Ai = Eu^r^ Hi = Vh & Ru- 
For j = 1, . . . ,m, define Wj = Bi . . . Bt such that Pj = O'' if i ^ Ci U . . . Cj, and Pj = Wu,ch H i = Ch G C„. 

Claim 4. E'i,...,E'^ and W[,...,W^ form a solution for M. 

Proof. Let I = li+. . .+lt. Let's consider a position in matrix M. Let P^- = vli . . . and Wj = Pi . . . P*. 
We discuss the following two cases. 

• Case 1. ^ P„ x (7„ for all u = 1, . . . , t. For each g <t, we have either Ag = O's or Bg = O'" since 
the submatrices Mi,...,Mt are independent. Therefore, E'^ AND Wj = O'. 
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• Case 2. (i, j) £ x C„. Since the submatrices Mi, . . . ,Mt are independent, we have (i, j) ^ Ry x 
for V u. We have AND By = O'" for each v ^ u.Thus, E'^ AND ^ O' if and only if 
Au AND Buy^O'-. 

□ 

By Claim 4, we have L{M) < L{Mi) + . . .+L{Mt). By Claim 3, we have L{M) = L{Mi) + . . .+L{Mt). □ 
Corollary G.4. There minimum length for the n x n diagonal-one matrix problem is exactly n. 
Theorem G.5. For most of n x n matrices, the minimum solution is at least ^ — n^^*^ for any fixed e > 0. 

Proof. Let e be an arbitrary positive constant. There are totally 2" many nxn binary matrices M. Assume 
TO < f — n^~'^- We consider how many matrices can be constructed by using /-list pairs with I < to. 

The total number of bits of a /-list pair is 2nl. The total number of /-list pairs is 2^"'. The total number 
of /-list pairs with / < to is J2"Li 2^"' < 2 • 2^"™ = o(2" ) since to < | — n^^^. Therefore, for most of n x n 
matrices, there is no ^ — n^"*^ solution. 

□ 

Theorem G.6. For the minimum length solution for a full rank matrix is greater than logn. 

Proof. Assume that there is a /-list pair Ei, . . . , E„i and Wi , . . . , Wn solution for an n x n matrix M of rank 
n. If there are two different i j with Ei — Ej, the i-th row and j-row are the same. Thus, the matrix is 
not full rank. If / = logn, then there is Ei to be all zeros. Thus the matrix M has a row to be all zeros. 
This makes M not to be full rank. Therefore, / > log n. 

□ 

Theorem G.7. For the minimum length solution for the nxn diagonal zero matrix is at greater than log n. 

Proof. The absolute value of the determinant of diagonal zero matrix is (n — 1). Thus, it is full rank matrix. 
It follows from Theorem IG. 61 □ 

Theorem G.8. Assume that k is an integer such that there is another integer 1 < /i < / such that (^'J > n. 
Then there is a /-list pair solution for the diagonal zero matrix problem. 

Proof. Let Ei, - ■ ■ ,i?„ be n different binary string of length / with exactly /i ones each. Let Wi be the 
complementary binary string of Ei. It is easy to see that Ei ■ Wi = and Ei ■ Wj > for i 7^ j. □ 

Corollary G.9. The minimum length for diagonal zero nxn matrix is between log n+l and log n + log log n. 

Proof. Assume that / is an even even number < logn -I- log logn. Let h = 1/2. By Stirling formula 
^ — „ , — = 1, we have (/ ) ^ . Thus, we can pick a / < log n + log log n such that (/ ) > n. The nxn 

V2irTi((^)") ' v'l' VttI — ° 00 vii/ — 

diagonal zero matrix is of rank n. The lower bound logn -I- 1 follows from Theorem |G.7[ 

□ 



G.2 Upper Bounds and Algorithm for Matrix Problem 

In this section, we show a randomized algorithm to handle a class of matrix problems. A matrix Mmxn 
has a list of independent submatrices Mi, • • • , Mt that contain all zero entries of M. We derive an upper 
bound of the solution for L{M) to be close to max(L(Mi), • • • ,L{Mt)). This implies the interesting bound 
max(L(Mi), • • • , L{Mt)) < L{M) < (1 + e) max(L(Afi), • • • , L{Mt)) + C)( '°g('^+") ) for an arbitrary positive 
constant e. 



We believe the following simple algorithm in Lemma G.IO for a random permutation is not new. For 
completeness, we include it here. 

RandomPermutation(n) 
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Let S* = {1,2, • • • ,n}; 
For i from 1 to n 

Select a random element from S 

Let 5* ^S-{a,}; 
Output aia2 ■ ■ ■ an', 
End of RandomPermutation 



Lemma G.IO. Every permutation of 1, 2, . . . , n has a equal probability to be generated by RandomPermu- 
tation(.) that runs in O(n^) time. 

Proof. Assume that aia2 • ■ ■ a„ be an arbitrary permutation. We just need to prove that Pi — aia2 ■ ■ ■ 0^0^+1 ■ ■ ■ a„ 
and P2 = aia2 ■ ■ ■ ai+ia^ • • • a„ have the equal chance to be generated. This is because a permutation can be 
converted into another permutation via a finite number of swaps between two neighbor items. Note that Pi 
can be converted into P2 by swapping the two elements and a^+i. Assume that the partial permutation 
aia2 ■ ■ ■ Ui^i have been generated by RandomPermutation(.). We have that 0^0^+1 and 0^+10; will be gener- 
ated by a equal probability RandomPermutation(.) to append to the last partial permutation 0102 • • • ai_i. 
The computational time follows from that fact that it takes 0(n) time to generate one element and update 
the set S in the algorithm. □ 

Lemma G.ll. Assume that 5*1 is a subset of ki elements of {1, 2, • • • ,n}. Let S2 be a subset of {1, 2, • • • ,n} 
and of size k2 with k2 > an for some a G (0, 1). Then with probability at most p < (1 — a)'^^ , P[5'2] n 5i = 
for a random permutation P of 1, 2, • • • , n, where P[S'2] is the set of the positions of elements ^2 in P. 

Proof. For two equal size subsets S and S" of {1,2,-- - the number of permutations P that make 



P[5'2] n = is the same as the number of permutations P that make P[5'2] n S" = 0. By Lemma G.IO 
the probability for P[52] n S* = is the same as the probability for P[S'2] n S" = 0, where P is a random 
permutation of 1, 2, • • • , n. Thus, we can assume that Si is {1, 2, • • • , ki}. 

The probability that none of the elements of ^2 are selected in the first fci positions is 



n n ~ 1 n ~ ki ~ 1 

Consider the function f{x) = ^'^^^^^ = We have its derivative function = -JT^^ < 0- 

Thus, f{x) is a decreasing function. We have < ^^^7^. Thus, 

n ^ k2 n — k2 — 1 n — fci — fc2 — 1 ,n ~ k2 ,1. 

, \ < ( (1) 

n n — 1 n — Ki — 1 n 

< (2) 

□ 

Lemma G.12. There is algorithm that given a binary matrix Mmxn, it gives a solution for the generalized 



matrix problem with S{M) length in 0{mn) times, where S{M) = min(m, n) as in Definition G.2 



Proof. Let M = {aij)mxn- Without loss of generality, assume n < m. For each column fc, let Wk — 
Qfe-i^Qii-fc^ For each row j, let Ej = Ej-.i . . . ej.„, where Cj^k = 1 if and only if aj^k = 1- D 

Definition G.13. For a strings, define a padding function pa(ii(s, fci, A;2, ^3) — O'"'^ I'^^sl'^^"'*' if |s| < fca, and 
padi(s, fci, fc2, fcs) is the empty string otherwise. We also define another padding function pad2(s, fci, fc2, fca) = 
lkiQk2 g^k3-\s\ j£ |g| <^ andpac?2(s, fci, fc2i ^3) is the empty string otherwise. For a permutationp = ii . . . i„ 
of 1, . . . , n, and a binary string s = ai . . . a„ of length n, define P{s,p) = ai-^ . . . ai^. 
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Lemma G.14. Let e be an arbitrary constant in (0, 1). Then there is an 0{{ni + n)^) time randomized 
algorithm such that given a binary matrix Mmxm a set H = {Mi, • • • , Mt} of independent submatrices that 
all zeros entries of M are in the submatrices of H, and also a solution with length li for each Mi £ H, it 
produces a solution for M with length (1 + e)U{H) + 0( '°g^"+'"^ ), where U{H) = max{/, : M, e H} is the 
largest length of the solution among all submatrices in H. 

Proof. Let fci = /j2 = ^U{H) + | log(m + n) and — U(H), where constant c is selected such that 

(1 _ I)! l°S("+m) < ^L. (3) 
4 Amn 

The total length for designing the matrix solution is z — ki + k2 + k^- We design Ei, - ■ ■ , and Wi , Wn 



of length z each for a solution for the matrix problem M. By Lemma G.12 we assume U{H) < min(m,n). 



For each submatrix in Af„ ~ {Ru, Gu) in H with Ru = {ii, ■ ■ ■ , is\ and C„ — {ji, . . . , jt}, let M„ have 
a solution e\'^\ . . . , eI"^^ and T4^j"' , . . . , wj"' , which has length = \e')^\ <U{H). Let p„ be a random 

permutation 1, 2, . . . , z. Now let Ei^ = P{padi{El^''),pu) for fc = 1, . . . , s and let Wj^ — P{pad2{Wj^^),Pu) 
for k — 1, . . . ,t. Thus, if row i and column j have G -Ru x C„ for some M„ = {Ru, Cy) G -ff, we 

have £;|"^ AiVi:» wj"^ = 0^ if and only if F(padi(£;f'^),p„) AiVDF(pad2(wj"^),p„) = 0^ if and only if 
Ei AND Wj = 0^ if and only if ^ = 0. 

For each row i, ii i ^ Rk for any Mk = {Rk,Ck) G H, then let Ei = V- . Similarly, for each column j, 
if i ^ Ck for any Mfe = {Rk, Cfe) G -ff, then let Wj — 1^. Thus, for row i and column j, if i ^ Rk for any 
Mk = {Rk,C'k) G i?, we always have AND Wj ^ 0^. Similarly, if j ^ Cfe for any Mk = (i?fc, Cfe) G H, we 
always have E, AND Wj ^ 0^ 



Let 



2fci + U{H) 

|C/(i?) + f log(m + n) 



2(f C/(iJ) + f log(m + n)) + U{H) 

eU{H) + ^ log(m 
2eU{H) + 2U{H) + f log(TO + n)) 
eC/(i7) + f log(m + n) 
4[/(7?) + ^log(m + n)) 
eU{H) + e • ^ log(TO + n) 
AU{H) + log(m + n)) 



(4) 
(5) 

(6) 



> .JJ. 4:._:L . (8) 

e 



(10) 



For an entry («,j), we consider the case that i G -Ru in Af„ = {Ru,Cu), and j G C^ in M^, = {Ry,Cy) 
for some u ^ v. For a sequence s = ai . . . a^, let Qi{s) be the the set of positions with bit 1 in s {Qi{s) = 
{i : tti — I and I < i < z}). Let -E^-"-*, • • • , E^"^^ and wj^\ • • • , wj"^ be a solution for submatrix -M„ with 
length lu = < U{H). Since i £ Ru, Ei = P{padi{E[^^),pu), where pu is a random permutation 

of 1, 2, • • • , z, and ife = i. Let • • • , e'^x] and ' ' ' j W^it''' be the solution for submatrix M„ with 

length ly = \e\"'^\ < U{H). Since j G Cy, Wj = P{pad2{Wy]^^),py), where py is a random permutation of 
1, 2, • • • , z, and yk — j- Since u ^ v, Pu and py are independent permutations. Qi{Ei) contains fci positions 
of {1,2,..., z}, and (5i(Wj) another fc2 positions of {1,2,..., z}. By Lemma G.ll with probability at most 



(1 - Qi{E,) n Qi(Wj) = 0. Therefore, with probability at most (1 - af\E, AND Wj = 
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Therefore, with probabiUty at most p — nm{l — a)^'^, there is a position («, j) in M with Ei AND Wj = 0^ 
not to be equivalent to a^j = 0. By the choice of c at equation ([s]), and inequahties Q to (10 1, we have 
small probability p defined above to be at most | . 

The computational time follows from the fact that for each Ei or Wj, we need at most 0((m + n)^) time, 
which is spent for generating a random permutation, from the solutions of matrices in H. □ 

Theorem G.15. Let e be an arbitrary constant in (0, 1). Assume that H is a set of independent submatriccs 
of M that all zero elements in M are in the submatrices of H, then V{H) < L{M) < (1 + e)V{H) + 
C>( '°g<"+"'^ ), where V{H) = max{i(M,) : M, e H}. 

Proof. It is trivial to see the inequality V{H) < L{M) since a solution for M automatically i mplies a solution 



for its submatrix. The inequality L{M) < (1 + e)V{H) + (9( '°s('^+"') ) follows from Lemma 



G.14 



□ 



Theorem G.16. Let e be an arbitrary constant in (0, 1). There is a randomized algorithm such that given 
a binary matrices Mmxn and a set H of independent submatrices of M, if all zero elements M are in the 
submatrices of H, then the algorithm returns a design with total length at most (1 + e)T(H) + )^ 
Furthermore, the time complexity of the algorithm is 0((n + m)^). 

Proof. It follows from Lemma |G.12| and Lemma |G.14| □ 

Lemma G.17. There is a 0(2^(^) (m + nY) time algorithm to find an optimal solution of length 

L{M) for a binary matrix M^xn- 



Proof. By Lemma G.12 it has a solution of length at most min(m, n). Assume m < n. We can find the least 
length by trying all numbers at most L{M). Each number takes at most 0(2^'*^)'") time. After vectors 
El, - ■ ■ , Em are fixed, it takes 0{mn) time to derive each Wi. This can be done to put the least number of 
zeros in the Wi to satisfy the zero entries of column i of Mmxn- Thus, it takes 0{mn?) time to derive all 
the columns after fixing rows El,-- - ,Em- Therefore, the total time is 0(2^(*^)(m + nf). □ 

Theorem G.18. Let e be an arbitrary constant in (0, 1). Then there is a randomized algorithm such that 
given a binary matrices M^xn and a set H of independent submatrices of M, if all zero elements M are 
in the submatrices of H , then the algorithm returns a design with total length at most (1 + e)V{H) + 
Qj- iog(»+™) where V{H) = max{i(Mi) : Mi G H}. Furthermore, the time complexity of the algorithm is 



Proof. It follows from Lemma G.17 and Lemma G.14 □ 



G.3 Revised Matrix Problem 

In this section, we revise the definition of the matrix problem, and show that how the optimal length of 
solution depends on the rank of the matrix. 

Revised Matrix Problem: Given a n x n nonnegative integer matrix M, find a list vectors Ei, . . . , E^ and 
a list vectors Wi , . . . , Wn of length I for each vector such that M{i, j) = Ei - Wj , where • is the inner product 
of two vectors. Minimize the length I. 

Theorem G.19. Let M be a n x n matrix with nonnegative elements. The minimum length solution for 
the revised matrix problem M is at least rank{M). 

Proof. Assume that there is a Mist pair Ei, . . . , En and Wi , . . . , Wn solution AI such that M(i, j) = Ei- Mj . 
Define Mt to be the matrix such that Mt{i,j) — Ei[t]Wj[j]. Mt is generated by the t-th bit of the Hist pair. 
We have M — Mi + M2 + . . . + M;. It is easy to see that rank of each Mt is at most one. 

By the well known fact in the linear algebra, we have rank(M) < rank(Mi)+rank(Af2)+. • .+rank(Af;) < /. 
Therefore, rank(Af) <l. □ 
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Corollary G.20. The minimum length solution for the full rank revised matrix with matrix size n x n 
problem is n. 

Corollary G.21. The minimum length solution for the revised matrix problem with matrix size n x n for 
diagonal one matrix is n. 

Corollary G.22. The minimum length solution for the revised matrix problem for diagonal zero n x n 
matrix is n. 

Corollary |G.22| is in contrast to Corollary |G.9[ 

Proof. The determinant of a diagonal zero matrix is n — 1. It is full rank matrix. It follows from Theo- 
rem [GH 

□ 
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